Non-Gaussianity and Dynamical Trapping in Locally Activated Random Walks 
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We propose a minimal model of locally-activated diffusion, in which the diffusion coefficient of a 
one-dimensional Brownian particle is modified in a prescribed way — either increased or decreased 
— upon each crossing of the origin. Such a local mobility decrease arises in the formation of 
atherosclerotic plaques due to diffusing macrophage cells accumulating lipid particles. We show that 
spatially localized mobility perturbations have remarkable consequences on diffusion at all scales, 
such as the emergence of a non-Gaussian multi-peaked probability distribution and a dynamical 
transition to an absorbing static state. In the context of atherosclerosis, this dynamical transition 
can be viewed as a minimal mechanism that causes macrophages to aggregate in lipid-enriched 
regions and thereby to the formation of atherosclerotic plaques. 
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Many-particle systems that consume energy for self- 
propulsion — active particle systems — have received 
growing attention in the last decade, both because of the 
new physical phenomena that they display and their wide 
range of applications. Examples include molecular mo- 
tors, cell assemblies, and even larger organisms [T]. The 
intrinsic out-of-equilibrium nature of these systems leads 
to remarkable effects such as non-Boltzmann distribu- 
tions P] , long-range order even in low spatial dimensions 
[3] and spontaneous flows jlj. 

At the single-particle level, the active forcing of a 
Brownian particle leads to non-trivial statistics. For ex- 
ample, it has been recently shown tS] [B] that a random 
walk which is reset to its starting point at a fixed rate has 
a non-equilibrium stationary state, as opposed to stan- 
dard Brownian motion. Another example is given by self- 
propelled Brownian particles [7] , which can yield sharply 
peaked probability densities for the particle velocity. 

In this letter, we consider a new class of problems in 
which the active forcing of a Brownian particle is localized 
in space. While the impact of localized perturbations on 
random walks has been investigated [8j, in part because 
of its relevance to a wide range of situations, such as 
localized sources and sinks [51 [TO]; trapping [TTJ [T2] or 
diffusion with forbidden [T3], hop-over [TJ] or defective 
[T5] sites, the role of local activation on Brownian-particle 
dynamics remains open. We present a minimal model 
of locally activated diffusion, in which the diffusivity of 
a Brownian particle is modified — either increased or 
decreased — in a prescribed way upon each crossing of 
the origin. 

A prototypical example is a bacterium in the presence 
of a localized patch of nutrients, which enhances the abil- 
ity of the bacterium to move, or, alternatively, toxins 
that impair bacterial mobility. This type of localized de- 
crease of mobility also underlies the dynamics of a cell 
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FIG. 1: a. Sketch of the different stages of atherosclerosis 
plaque formation: (1) rapid diffusion of a "free" macrophage 
cell; (2) upon entering a localized lipid-enriched region, the 
macrophage accumulates lipids, and thereby grows and be- 
comes less mobile; (3) after many crossings of the lipid- 
enriched region, the macrophage eventually gets trapped, re- 
sulting in the formation of an atherosclerotic plaque, b. 
Sketch of a one-dimensional particle trajectory of the model 
of locally decelerated random walk. 



(e.g., a macrophage) that grows by accumulating smaller 
and spatially localized particles, such as lipids (Fig. [T]) 
[Tni[T7]. As the cell grows, its ability to move decreases 
and the ultimate result is the formation of an atheroscle- 
rotic plaque [18| . The spatial localization arises from 
the presence of lipids at specific points in the arterial 
network; these lipids can be located, as is now well ac- 
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cepted, by the properties of the blood flow (TH] • Obser- 
vations show that macrophages that have accumulated 
lipids move more slowly. Eventually the macrophage 
stops in an lipid-enriched region, resulting in the for- 
mation of an atherosclerotic plaque [201 HI]- Here we 
propose a simple model to account for this local mobil- 
ity decrease and address the particular questions of (i) 
the potential trapping of cells in locally lipid-enriched 
regions, and (ii) the kinetics of the resulting segregation 
process when it exists. 

Our formalism allows us to describe both the situations 
of decreased and increased localized mobility changes. 
We show that this type of perturbation has remarkable 
consequences on the diffusion process at all scales. We 
stress that the diffusion coefficient of the active particle 
at any time depends on the entire history of the tra- 
jectory. Thus the evolution of the particle position is 
intrinsically non-Markovian [22'-'25' . Our main findings 
are: (i) The probability distribution of the position has a 
non-Gaussian tail, (ii) For local acceleration, a diffusing 
particle is repelled from the origin, so that the maximum 
in the probability distribution is at non-zero displace- 
ment, (iii) For local deceleration, a dynamical transition 
to an absorbing state occurs. For sufficiently strong de- 
celeration, the particle can get trapped at the origin at a 
finite time. The exact time dependence for the particle 
survival probability is determined explicitly. Conversely, 
if the deceleration process is sufficiently weak, the parti- 
cle never gets trapped. This dynamical transition to an 
absorbing state provides a minimal mechanism that could 
help understand the formation kinetics of atherosclerotic 
plaques. 

The Model. A one-dimensional diffusing particle is ac- 
celerated or decelerated whenever it crosses the origin 
a; = according to the following Langevin equations: 



D = fiD)S{x), 



(1) 



where ^ is a Gaussian white noise of intensity one, D the 
particle diffusion coefficient, x the particle position, S{x) 
the Dirac distribution, and f{D) an arbitrary prescribed 
function that accounts for the local activation. For sim- 
plicity, we assume that the particle is initially at x = 
with D — Dq > 0. Note that: (i) both the position x and 
the diffusion coefficient D are random variables; (ii) as 
mentioned previously, the evolutions of a; or _D alone are 
non-Markovian; (iii) the function f{D) can be positive 
(local acceleration) or negative (local deceleration), but 
with /(O) = so that D remains non negative. 

Following standard steps, the corresponding Fokker- 
Planck equation for the joint distribution of position 
X and diffusion coefficient D at time t, P{x, D.,t), is: 



dP _ d^P 
~dt ^ 



Si^)^^^^^ - mSix)SiD), (2) 



where the last term of the right side accounts for the ab- 
sorbing state at {x = 0, D = 0). The explicit expression 
for X{t) is determined demanding that P is normalized. 



from which we obtain 



A(t)= hm [/(i?)P(0,D,t) 



(3) 



When f{D) is positive, then D is always non zero. In 
this case, the particle is never trapped and A(t) = at 
all times. While intuitively obvious for local acceleration 
{f{D) > 0), we show below that X{t) can equal zero for 
local deceleration processes. 

Local Acceleration, f{D) > 0. Laplace transforming 
Eq. ([2| gives 



sP + D 
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where P = P{x, D, s) is the Laplace transform of the 
probability distribution. For x ^ 0, the solution is 



P{x,D,s) = A{D,s)e- 



(5) 



where the coefficient A{D, s) is determined by integrat- 
ing Eq. ^ across a; = to obtain the jump of the first 
derivative of P with respect to x at this point: 
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dP 
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dP 

x=o+ dx 



d[fP{x = 0)] 
dD 



-S{D-Do). 



Using Eq. ([5|, we have 
.dA 



f 



dD 



A = S{D-Do) 



(6) 



When f{D) is positive, then A{D,s) = for D < Dq, 
while for D > Dq the solution to ([6| is 



(7) 



with 



F{D) 



D 



'D' 



D„ \f{D')\ 



dD'. 



The unknown function B{s) is determined by the jump 
of A at Dq: 



AiD+,.s)~A{DQ,s) = 
which finally yields 

P{x,D,s)^eiD-Do) ^ 



1 



f{D) 



fiDo) ' 



(8) 



where O is the Heaviside step function and 

-2F{D). 



Laplace inverting this expression, we obtain the joint dis- 
tribution 



P{x,D,t) = Q{D-Do) 



f{D)V^ 



(9) 



3 



The marginal distribution with respect to x, that is, 
the probability distribution of positions, is obtained by 
integrating Eq. ^ over all D in the range [Dq, oo]. While 
it does not seem possible to evaluate this integral ana- 
lytically, the large- a; behavior can be obtained by the 
Laplace method. For the illustrative case where f{D) is 
a constant (that we define as a), this method gives 



exp 



8|x|3/2 



oo , 



(10) 



which we numerically checked is close to the exact value 
P{x,t) over a wide spatial range. We wish to empha- 
size two important features of this result for P{x,t) that 
are in marked contrast with the Gaussian propagator of 
the usual Brownian motion: (i) P{x, t) generally has a 
non-Gaussian tail; (ii) P{x,t) reaches its maximum at 
a non-zero displacement. Equation (10) shows that the 
location of this maximum asymptotica ly grows as t^^^ 
when f{D) = a. Thus local acceleration pushes a diffus- 
ing particle away from the origin. 

From the general expression Q, the marginal distri- 
bution with respect to D can also be easily obtained by 
integration over x. We find 



P{D,t)^e{D^Do) 



2VD 



fiDWnt 



(11) 



In the particular case of /(D) — a, Eq. (11) shows that 



the diffusion coefficient of the particle asymptotically 
grows as i^/^. 

As a byproduct, Eq. ( |Tl| ) also provides the distribution 
of the local time T{t) spent by the particle in the active 
zone (the origin for the present case) up to time t. Using 
the second of Eqs. ([T]), this basic observable in the the- 
ory of diffusion, which has dimension of time per unit of 
length [37], is related to the diffusion coefficient at time 
t by 



-it) 



* r° dD' 



(12) 



Thus the distribution of the local time, defined as 7^(r, i), 
is given by V{T,t) = f{D)P{D,t), with P{D,t) given 
by Eq. (11) and D implicitly defined as a function of r 
in Eq. (12). For the illustrative case of f{D) — a, the 



distribution of the local time at time t therefore is 

.3/2] 2' 
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exp< 
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.(13) 

This result strongly contrasts with the Gaussian distri- 
bution that arises in the case of Brownian motion, which 



can be recovered from Eq. ( 13 1 in the limit a — >■ 0: 



(14) 



Notice, in particular, the typical local time for an accel- 
erated particle with f{D) — a grows as t^^^ instead of 
i^/^ in the case of Brownian motion. 



It is worth noting an intriguing dichotomy with a 
discrete-time version of local acceleration — the "greedy" 
random walk [25] • In this discrete model, the step length 
£k after the k*"^ return of a random walk to the origin is 
given by £k — k". To match with the continuous model 
with f{D) = D", one must choose the value a = 1/2. 
With this choice, Eq. (13) of [55] gives, ignoring all multi- 
plicative factors, P{x,t) oc cc^/^/i exp[— x^'^'^/i], which is 
different from ( 10 ) . The source of this discrepancy is that 



the probability of being at the origin is not affected by 
the enhancement mechanism of greedy walks [28] , while 
this return probability is fundamentally modified in the 
case of locally-activated random walks, as seen explicitly 
from the distribution of the local time (13 1. Thus our 



locally activated diffusion model cannot be viewed as the 
continuous limit of the greedy random walk. However, it 
can be shown that it is the continuous limit of a discrete 
space and continuous time random walk whose jump fre- 
quency is modified at each visit of the active site, which 
is intrinsically different from the greedy random walk. 

Local Deceleration, f{D) <0. Following the same anal- 
ysis as that used for local acceleration, the Laplace trans- 
form of the joint distribution is 



(15) 



where A is the Laplace transform of A(<) defined in 
Eqs. ([2| and ([3|. Using these defining relations for X{t), 
Eq. (fl5|) gives 



X{s)= ]im[f{D)P{0,D,s)] 

D^O 



/4s F 



(16) 



where we define 



F{D)^ 



'D' 



IfiD') 



dD' . 



In this result for A, we have used 6{D)f{D) — 0, since 
/(O) = by the definition of our model. The important 



feature of Eq. ( 16 1 is that A(s) = as soon as F diverges. 
Thus our final result is 



Pix D ,^ Q(^o-^) ,^z./i , 5{x)5{D) ^^^^ 



which gives, after Laplace inversion. 



(17) 



Pix,D,t): 



\f{D)\ 



T{t)S{x)S{D). 



Here 



/ 1 \ID' 



(18) 



(19) 



is the trapping probability, namely, the probability that 
the particle becomes stuck at a; = by time t because 
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the diffusion coefficient has reached zero. As a corol- 
lary, the survival probability is given by S{t) = f — T(<), 
and we have obtained this quantity for an explicitly non- 
Mar kovian process. We also mention that, as in the case 
of local acceleration, the joint distribution easily gives the 
marginal distributions of the position and the diffusion 
coefficient, as well as the local time. 

A fundamental consequence of the local deceleration 
of a Brownian particle is that two different dynamical 
regimes emerge. We illustrate these regimes for the par- 
ticular case where f{D) = -D^ as D ^ 0. If the decel- 
eration is sufficiently strong, which occurs when /3 < 3/2, 
there is a non-zero probability for the particle to get 
trapped at the origin. More precisely, the survival prob- 
ability has the asymptotic behavior 



S{t) 



3/2- 



7rt(3- 2/3) 



(20) 



Thus in this regime of strong deceleration, the survival 
probability has the same scaling with time as in the case 
of a usual Brownian particle in the presence of a perfect 
trap. In the opposite case of B > 3/2, then S'(f) = 1 
for all f > and the particle never gets trapped at the 
origin. Thus a locally decelerated Brownian particle un- 
dergoes a dynamical transition to the absorbing state 
{x = Q,D = Q) as the deceleration strength increases. 
Mathematically, this transition occurs at the point where 
F is no longer divergent. 

In conclusion, we introduced a minimal model of lo- 
cally activated diffusion, in which the diffusion coeffi- 
cient of a Brownian particle is modified in a prescribed 
way at each crossing of the origin. In one dimension, a 
purely diffusing particle hits the origin of the order of ^/t 
times after a time t. Consequently, the local activation 



mechanism is repeatedly invoked during the trajectory 
of a Brownian particle. Thus the asymptotic dynam- 
ics of a Brownian particle is globally affected, leading to 
markedly different behavior than that of pure diffusion. 
Since the unusual properties of local activation rely on 
the recurrence of Brownian motion, we anticipate that 
qualitatively similar, but quantitatively distinct, behav- 
ior would arise in two dimensions. 

Our model encompasses both the situations where the 
Brownian particle is locally accelerated or decelerated. 
For local acceleration, the probability distribution is non- 
Gaussian and multi-peaked, with maxima away from the 
origin no matter how weak the acceleration. For suffi- 
ciently weak local deceleration, a Brownian particle man- 
ages to avoid getting trapped at the origin in spite of its 
recurrence. However, for strong deceleration, there is a 
dynamical transition to an absorbing state in which the 
particle ultimately gets trapped at the origin. 

In the context of atherosclerosis mentioned initially, 
the dynamical transition to an absorbing state can be 
viewed as a minimal mechanism that leads to the segre- 
gation of macrophages in lipid-enriched regions and thus 
to the formation of atherosclerotic plaques. Our model 
suggests that even in absence of chemical signals (such 
as chemokines or cytokines) that can bias the motion of 
cells, there exists a critical intensity of the mobility de- 
crease, which depends on the local lipid concentration, 
beyond which an atherosclerotic plaque will occur. Our 
model can also help understand the kinetics of this plaque 
formation. 
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